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FOREWORD 

The  sub-project  was  funded  by  a Department  of  Defense  Themis 
Contract  No.  NOOOll*-C8-A-Ol60 , which  was  under  the  administration  of 
the  Institute  of  Atmospheric  Sciences  and  under  the  suoervision  of 
Dr.  R.  C.  Huntsingcr.  Contributors  were  Mr.  T.  E.  Warborg,  Mr.  G.  K. 

Bien  and  Mr.  T.  P.  Nicholas,  South  Dakota  School  of  Mines  and  Technology. 
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ABSTRACT 
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Over  a four  year  period,  from  September  10oT  to  July  1971,  the  rate 
of  growth  of  a 3Qu  radius  drop  to  50  or  lOQy  radius  by  a process  known 
as  accretion  has  been  studied  as  part  of  the  numerical  modeling  effort 

of  Project  Themis.  A mathematical  model  was  derived  which  considers  a 
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single  3Qp  radius  drop  falling  at  its  terminal  velocity  through  a field 
of  smaller,  5 to  20y  radius  droplets,  and  growing  by  colliding  and 
coalescing  with  droplets  in  its  path.  The  droplets  also  fall  at  their 
terminal  velocities  and  the  rate  of  growth  of  the  dron  is  proportional 
to  the  difference  in  the  terminal  velocities  of  the  colliding  drop  and 
droplet  as  well  as  the  area  of  sweep  of  the  falling  drop  and  the  number 
of  collectible  droplets  per  unit  volume  of  cloud.  Hydrodynamic  and 
electrostatic  forces  influence  the  rate  of  growth  and  are  accounted 
for  by  introduction  of  a '^collision  efficiency!'  term  to  the  analysis. 
Electrical  forces  also  influence  the  terminal  settling  velocities  of 
the  particles,  *so  Etoke's  law  of  drag  on  a sphere  and  the  intermediate 
law  were  derived  incoroorating  these  forces.  . 

V7arborg  (l)  studied  the  effect  on  drop  growth  rate  of  0 to  3600 
V/cm  electrical  field  strengths,  10-18  to  10“* 3 coulomb  charges  on  the 
drop  and  droolets  and  various  populations  of  droplet  sizes  through 
which  the  drop  falls.  A 5u  radius  monodisoersion  of  droplets  was 
considered  as  well  as  a pair  of  natural  cloud  droplet  distributions. 

His  study  emphasizes  the  imoortance  of  the  natural  droplet  distribution 
and  the  influence  of  electrical  effects  on  the  system  as  a 30y  radius 
drop  grows  to  50u  radius . 
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Bien  (3)  extended  the  growth  range  of  the  dron  to  100u  radius 
using  the  unpublished  collision  efficiencies  of  fartor.  He  investi- 
gated the  effect  on  dron  growth  rate  of  field  strengths  of  298  to 
2Q80  V/cm  and  droolet  charges  of  0 to  10r2  statcoulorabs , assuming 
that  the  charge  on  the  droplet  is  at  all  times  proportional  to  its 
surface  area.  He  studied  a number  of  monodispersions  and  natural 
droplet  distributions  and  approximated  the  droolet  distribution 
with  a mean-sized  ronodisoersion . 

Nicholas  (5)  considered  the  case  of  allowing  the  drop  to  collide 
with  droplets  ranging  from  1 to  lOOp  radius  some  of  which  could  be 
larger  than  as  well  as  smaller  than  the  drop.  This  reduces  the  time 
necessary  for  the  drop  to  grow  to  100p  radius  by  as  much  as  78? . 
Finally,  the  droplet  distribution  was  allowed  to  evolve  in  a theoreti- 
cal manner  toward  a situation  in  which  each  size  category  would  have 
the  sane  number  of  droplets  per  unit  volume  of  cloud.  The  effect  was 
again  a significant  reduction  in  necessary  growth  times. 

The  results  of  all  workers  indicate  that  the  effect  of  an  imposed 
electric  field  and  charges  on  drops  significantly  increase  the  growth 
rate  of  the  dron  once  minimum  values  of  these  electrical  forces  are 
reached. 

In  all  cases  studied,  increasing  the  electrical  influences  in 
the  model  resulted  in  a decrease  in  the  time  for  the  30p  drop  to  grow 


to  100u  radius. 
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Bien  found  that  varying  the  liquid  water  content  of  the  droplet 
distribution  has  a direct  effect  on  the  growth  rate  of  the  drop.  He 
indicated  that  the  drop  growth  rate  is  a linear  function  of  the 
available  water  in  the  cloud. 

Research  efforts  concerning  modification  of  convective  clouds 
have  resulted  in  a mathematical  model  of  cloud  drop  growth  which  can 
closely  approximate  conditions  existing  in  an  actual  cloud.  'todifi- 
cations  have  been  made  to  the  model  with  the  goal  of  improving  the 
reliability  of  the  evaluations  of  cloud  drop  growth  by  accretion 
under  various  electrical  influences . 
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SUMMARY  ADD  CONCLUSIONS 


One  of  the  major  processes  by  which  cloud  water  in  transformed 
into  rainwater  is  known  as  accretion.  One  type  of  accretion  results 
from  collisions  between  droplets  subjected  to  external  forces  such  as 
gravity  and  electric  force  fields.  The  probability  of  two  droplets 
coalescing,  combining  upon  collision,  was  the  subject  of  the  research 
of  Sub-project  III  of  Project  Themis  at  the  South  Dakota  School  of 
Mines  and  Technology.  A mathematical  model  has  been  developed  investi- 
gating a 30p  radius  drop  falling  at  its  terminal  velocity  through  a 
field  of  smaller  droplets  (5  to  20u  radius),  also  falling  at  their 
respective  terminal  velocities.  The  drop  grovs  by  overtaking  droplets 
in  its  path,  colliding  with  and  possibly  combining  with  them.  Warborg 
(l)  has  derived  the  growth  equation  which  is  based  on  a simple  con- 
servation of  mass  principle.  The  rate  of  change  of  drop  radius. 


dr 

dt 


-JT_ 

3r2 


rl=5 


(1) 


The  summation  accumulates  the  individual  contributions  of  each 
droplet  size  category,  rj , to  the  increase  in  dron  radius,  r.  The 
remainder  of  the  factors  are  explained  in  the  following  sections. 

In  a vacuum,  the  maximum  cross-sectional  area  of  sweep,  A Cm2, 
for  which  the  drop  which  would  Just  contact  a droplet  would  be 
ir(r  ♦ x*j ) 2 , where  r and  rj  are  the  radii  of  the  dron  and  droplet  in 
centimeters.  Upon  introducing  hydrodynamic  forces  and  other  external 
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electrical  influences  such  as  charged  particles  in  an  electrical  field, 
the  theoretical  cross-section  of  sweep  can  be  considered  to  be  modified 
by  a correction  factor  termed  the  "collision  efficiency". 

!fodification  of  equation  (l)  to  accommodate  the  collision 
efficiency  term  yields: 


dr 

dt 


20 


z 

ri=5 


r^n  Ec  (r  + rj  )2  Ut 


(2) 


Warborg  used  the  collision  efficiencies  calculated  by  Plunlee  which 
provided  values  for  collision  between  drons  less  than  50v  radius  under 
various  electrical  influences  while  Bien  (3)  and  Nicholas  (5)  used  the 
theoretically  determined  collision  efficiencies  of  Sartor  in  their 
analyses . 

Computation  of  terminal  velocity,  cm/sec,  for  falling  drops  has 
been  obtained  by  modifying  the  Stokes'  equation  for  a suhere  falling 
through  a viscous  medium,  to  include  a term  accounting  for  electrostatic 
force  on  a char/red  dron  in  the  presence  of  an  electric  field.  Bien  (3) 
shows  this  derivation  in  his  master's  thesis  which  holds  for  Reynolds 
numbers  less  than  2.  The  final  enuation  becomes: 


IN 


2r2 

9u 


(p  - p) 
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3qE 
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(3) 


For  Reynolds  numbers  between  2 and  500,  another  form  of  the  equation 
from  McCabe  and  Smith,  see  Nicholas  (5),  called  the  intermediate  law. 


is  used.  The  formulation  is: 


For  the  distribution  of  droplets,  N through  which  the  drop  falls , 
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Warborg  (l)  used  a lOp  monodisoersion  and  a natural  disnersion  of  5 to 


20p  droplets  determined  for  fair  weather  cumulus  clouds  by  Battan  and 
Reitan  shown  in  Figure  (l).  He  indicated  that  the  full  droplet  dis- 
tribution significantly  increases  the  growth  rate  of  the  dron  as 
compared  with  the  monodispersion.  Bien  (3)  used  the  same  natural  droplet 
distribution  and  also  investigated  various  other  monodispersions.  His 
analysis  shows  that  a dronlet  distribution  can  be  approximated  by  a 
mean-size  monodispersion  with  the  same  liauid  water  content.  A.  mean 
droplet  radius  was  used  to  determine  the  size  of  the  representative 
monodispersion.  The  equation  takes  the  form: 


Where  rx  is  the  mean  radius  of  the  monodispersion  in  cm,  representative 
of  the  full  dronlet  distribution.  The  effect  of  aoolying  this  principle 
to  the  model  is  illustrated  in  Figure  (3).  Nicholas  (5)  then  extended 
the  size  range  of  the  droplet  distribution  by  introducing  a new  distri- 
bution which  was  presented  by  Junge  including  1 to  lOOp  radius  droplets. 
The  distribution  is  shown  as  curve  tl  in  Figure  (2).  Comparison  of  the 
growth  rate  produced  by  this  distribution  with  that  produced  by  the 
Battan  and  Reitan  distribution  shows  the  effect  of  allowing  the  drop  to 
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oW  tv  collecting  droplets  which  may  be  larmer  as  well  as  smaller  than 
itself* 

Comparison  of  the  constant  spectrum  curves  of  Figures  (5)  and 
(7)  and  Figures  (6)  and  (8)  illustrates  that  using  the  Junge  distribution 
results  in  substantially  reduced  growth  times  as  compared  to  the  cases 
f involving  the  Battan  and  Reitan  distribution.  At  small  electrical 

parameter  values  the  percent  of  decrease  in  the  growth  times  reached 
as  high  as  20Q1?.  At  larger  electrical  influences,  this  difference 
becomes  less  oronounced  reaching  a minimum  of  approximately  100!?  at  the 
maximum  values  of  the  field  strength  and  surface  charge  density  studied. 

Finally,  Nicholas  (5)  developed  a theoretical  time  evolution  of 
the  droplet  spectrum  based  on  the  limit  of  equal  number  density  for 
each  droplet  size  at  some  future  time.  The  time  evolutions  for  the 
Battan  and  Reitan  and  the  Junge  distributions  are  depicted  by  Nicholas 
(5)  as  Figures  (l)  and  (2).  The  original  distribution,  indicated  by 
tl  is  allowed  to  evolve  toward  a limiting  condition  represented  by 
tire-curve  tU  wherein  all  of  the  droplet  size  categories  contain  an 
ccual  number  of  droplets.  Intermediate  values,  curves  t2  and  t3  are 
obtained  by  linear  interpolation  between  points  on  tl  and  tU  for 
constant  droplet  radius.  Comparison  of  "constant  spectrum"  and 
evolvinr  spectrum"  curves  in  Figures  (5)  to  (8)  indicate  that 
rubstantial  growth  time  savings  (22  to  78!?)  , result  from  allowing 
-f*e  distributions  to  age  with  time.  Inconsistent  portions  of  Figures 
.”)  and  (8)  at  low  surface  charge  density  resulting  in  a lower  value 
'-r  rro*..-th  time  of  the  constant  spectrum  have  been  traced  to  the 
"-lision  efficiency  data  of  Sartor.  Incorporation  of  the  evolving 
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droplet  spectrum  into  the  model  produced  a si/mificant  increase  in 
the  growth  rate  of  the  drop  for  almost  every  case  tested. 

Electrical  effects  have  been  central  to  all  three  analyses  and 
have  been  used  throughout  as  parameters  of  the  commuter  simulations 
of  the  mathematical  models.  Warborg,  Bien,  and  Nicholas  all  depicted 
growth  times  for  a specific  electrical  field  strenirth  as  a function 
of  the  charge  of  the  droplet  of  the  system.  The  assumption  is  made 
that  the  charge  on  a drop  is  at  all  times  proportional  to  the  surface 
area  of  the  drop,  or  the  surface  charge  density  is  constant  for  one 
computer  run.  Hence,  the  value  in  the  abscissa  of  Figures  (5)  through 
(8)  for  instance,  multiplied  by  Uirr2  gives  the  charge  in  st  at  coulombs 
on  the  drop  or  radius  r.  The  constant  electric  field  strenirth  is 
marked  on  each  graph.  In  all  cases  observed,  an  increase  in  the 
electric  field  strength  produced  a corresponding  increase  in  the 
growth  rate  of  the  drop.  Likewise,  increasing  the  surface  charge 
density  on  the  dronlets  increases  the  drop  growth  rate  correspondingly. 
This  trend  agrees  with  other  work  in  the  field  and  substantiates  the 
assumptions  made  in  the  mathematical  model. 

Bien's  (3)  analysis  determined  a linear  dependence  of  the  drop 
growth  rate  on  cloud  liquid  water  content.  The  dependence  is  illustrated 
in  Figure  (U). 

The  direction  of  the  research  efforts  regarding  the  mathematical 
model  has  been  toward  formulation  of  a more  realistic  drop  growth  model, 
including  effects  observed  in  an  actual  cloud  parcel.  Knowledge  of  the 
growth  characteristics  of  cloud  droplets  and  the  factors  which  influence 
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the  growth  c.  .1  provide  valuable  insight  into  the  understanding  of  the 
overall  precipitation  mechanism  in  the  modification  of  convective 
clouds . 


11 


} 


TIO,'r‘:!CLA"’TTT’r 


A 

E 


E 

c 


R 


P,  Pp 


Cross-sectional  area  of  sweep  of  a fall  in  r drop,  cm*1 . 
Electric  field  strength,  V/cm. 

Collision  efficiency,  dimensionless. 

Acceleration  of  aravitv,  9^1  cr./sec2. 

Number  of  dronlets  of  radius  r per  cm3  of  cloud. 
Droplet  charge,  st at coulomb . 

Dron,  droplet  radii,  cm. 

Tine,  sec. 

Approach  velocity  of  large  drop  to  small  droplet, 

( TJ  - U ) cm/sec. 

Li  s 

Dynamic  viscosity  of  air,  g/cn/sec. 

Densities  of  air  and  water  respectively,  g/cm3 . 
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Figure  (l). 


Figure  (2). 


Figure  (3). 


Figure  (U). 


Figure  ( 5 ) . 


Figure  (6). 


Battan  and  Reitan  discrete  droplet  distribution  and 
theoretical  evolution,  tl  = original  distribution, 
t2  -*■  tl*  = evolution  to  constant  distribution  function. 
Liquid  water  content  = 0.665  g/m3. 

Junge  droplet  distribution  end  theoretical  evolution, 
tl  = original  distribution,  t2  ■+  tU  = evolution  to 
final  distribution.  Liquid  water  content  = 0.U21  g/m3 . 
Growth  times  for  drop  to  grow  from  30  to  100  y radius 
versus  drop  charge  for  5,  10,  and  18.32  y radius  mono- 
dispersions  for  the  Battan  and  Reitan  distribution. 

Curve  marked  "28.32  y Jtonodispersed  Cloud"  represents 
calculated  mean  monodispersion.  Field  strength  E = 

298  volts/cm. 

Growth  times  for  drop  to  grow  from  30  to  100  y radius 
versus  drop  charge  for  5 , 10 , 15 , and  20  y radius 
monodispersions  for  the  Battan  and  Reitan  droplet 
distribution.  Field  strength  E = 2980  volts /cm. 

Growth  times  for  drop  to  grow  from  30  to  100  y radius 
versus  drop  charge  for  constant  and  evolving  Battan  and 
Reitan  distribution.  Field  strength  E = 298  volts/cm. 
Growth  times  for  drop  to  grow  from  30  to  100  y radius 
versus  drop  charge  for  constant  and  evolving  Battan 
and  Reitan  distribution.  Field  strength  E = 2980  volts/cm. 
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Figure  (7). 


Drop  growth  times  versus  drop  charge  for  constant 
and  evolving  Junge  distribution.  Field  strength  E 
298  volts /cm. 

Figure  (8).  Drop  growth  times  versus  drop  charge  for  constant 
and  evolving  Junge  distribution.  Field  strength  E 
2930  volts /cm. 
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